Phase-Dependent Electronic Specific Heat in Mesoscopic Josephson Junctions 
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We study the influence of superconducting correlations on the electronic specific heat in a diffusive 
superconductor-normal metal-superconductor Josephson junction. We present a description of this 
system in the framework of the diffusive-limit Green's function theory, taking into account finite 
temperatures, phase difference as well as junction parameters. We find that proximity effect may 
lead to a substantial deviation of the specific heat as compared to that in the normal state, and that 
it can be largely tuned in magnitude by changing the phase difference between the superconductors. 
A measurement setup to confirm these predictions is also suggested. 

PACS numbers: 74.78.Na,74.25.Bt, 74.45. +c,73.23.-b 



Superconducting circuits attract much interest thanks 
to their potential in applications such as single-electron 
devices [ll, I3l and Josephson junctions [3], as well as for 
radiation detection Q. In this context, proximity effect 
in mesoscopic normal metal-superconductor nanos- 
tructures has drawn researchers attention @, 0| ■ In such 
systems the normal metal can acquire superconducting- 
like properties which manifest themselves in macroscopic 
observables, such as the Josephson current Q- Although 
the focus has been mainly given to the electronic proper- 
ties of proximity structures, thermal transport is recently 
under the spotlight fsl. Yet, very little is known about 
the thermodynamic quantities of these systems, for in- 
stance, the entropy as well as the specific heat, which 
could be crucial for future applications in nanoelectron- 
ics. 

In this Letter we theoretically address the influence 
of superconducting correlations on the electronic spe- 
cific heat of a diffusive superconductor-normal metal- 
superconductor (SNS) Josephson junction. We show 
that, thanks to proximity effect, the specific heat of the 
N region can be dramatically different from that in the 
absence of superconductivity. In particular, at low tem- 
perature, it is substantially suppressed with respect to 
the normal state, and turns out to be largely tunable in 
magnitude by changing the phase difference between the 
S reservoirs. Such peculiarity of mesoscopic SNS Joseph- 
son junctions may have impact for the implementation of 
novel devices, for instance, ultrasensitive single-photon 
detectors based on proximity effect P . 

The mesoscopic system under investigation [see the in- 
set of Fig. [T]Ja)] consists of a diffusive normal metal 
wire of length L in contact with two superconducting 
electrodes kept at zero potential, thus defining a SNS 
Josephson junction. The wire is oriented along the x 
direction, and we assume its transverse dimensions to 
be much smaller than L, so that the wire can be con- 
sidered as quasi-one-dimensional. The contact with su- 
perconducting leads induces in the N region correlations 
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FIG. 1: (color online) (a) Comparison of the SNS electronic 
specific heat C vs temperature T with that in the normal state 
(N) and in the superconductor (S) calculated for Aq = 10i?Th 
and (j) = 0. The inset shows a scheme of the SNS junction of 
length L. X — denotes the middle of the junction, and the 
system is assumed quasi-one-dimensional, (b) c vs T at dif- 
ferent positions x along the N wire calculated for Aq — lO-Exh 
and = 0. Also shown is the specific heat in the absence of 
superconductivity (N). Tc is the critical temperature of the 
superconductor. 

through proximity effect Q, which is responsible for the 
Josephson supercurrent flow across the structure. 

The proximity effect in the normal metal region can 
be described with the Usadel equations [13] ■ According 
to the parametrization of the Green's functions given in 



Ref. the Usadel equations in the N region can be 



written as [10|, [11 1 



hDdle = ~2iEsmh{0) + ^ (d^xf sinh(26l) 

sinh{2e)d^0d,X + smhHO)dlx = 0, (1) 

where D is the diffusion coefficient and the energy E is 
relative to the chemical potential in the superconductors. 
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9 and x are, in general, complex scalar functions of posi- 
tion (x) and energy. For perfectly transmitting interfaces 
fl^ the boundary conditions at the NS contacts reduce 
to 



e{±L/2) 

x{±L/2) 



arctanh(A/£') , 
±(f>/2, 



(2) 



where is the phase difference between S electrodes, and 
A is the superconducting order parameter. We adopt a 
step- function form for the order parameter, i.e., constant 
in S and zero in the N wire, although A is in princi- 
ple x-dependent and can be determined self-consistently 
Furthermore, we assume the usual BCS tempera- 
ture dependence of A(T) with critical temperature Tc — 
Ao/(1.764fcB), where Aq is the zero-temperature order 
parameter, and /cb is the Boltzmann constant. The den- 
sity of states (DOS) in the N region is then given by 

JV{x, E, T, 0) = TVpRe {cosh [0{x, E, T, 0)] } , (3) 

where A/p is the normal-state DOS at the Fermi level. 
Assuming that the S electrodes are not influenced by the 
proximity effect, the total electronic specific heat (C*°') 
of the SNS junction is determined by 



c*°'(r. 



C(T,0)-f Cs(T), 



(4) 



where Cs(T) is the temperature-dependent specific heat 
of the superconducting leads, while 



dS{T, , 
dT 



(5) 



is the specific heat of the N region. Furthermore, the 
entropy S{T, cj)) is defined as [14 1 



5(r,, 



dxdEM{x,E,T,(j)) (6) 



4 i-L/-^ 

L J-L/2 Jo 

X [f{E) Mf{E)) + (1 - f{E)) ln(l - fiE))] , 



and f{E) = (1 + exp[E/{kBT)])-^ is the Fermi-Dirac 
quasiparticle distribution function at temperature T. 
The characteristic energy scale for the N region is ETh — 
HD/L"^, usually referred to as the Thouless energy. In 
the following we shall focus on long SNS junctions only 
(Aq ^ -ETh), as it is the relevant regime for metallic 
diffusive Josephson junctions [lilllit. 

By numerically solving Eqs. ([1]) with the above given 
boundary conditions [see Eqs. ^] we can determine the 
DOS as a function of energy and position for a given (j) 
and temperature. In particular, the DOS is an even func- 
tion of energy, and exhibits a minigap (Eg) for \E\ < Eg 
17| which depends, in general, on i^xh and (p. In addi- 
tion. Eg w 3.2£^Th for long junctions ai (f> = [l^, and 
decreases by increasing 0, vanishing at = tt. Moreover, 
the density of states at NS interfaces (i.e., at x = ±L/2) 
equals the BCS superconducting DOS, while it tends to 
Afp for \E\ > Aq. 
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FIG. 2: (color online) (a) Specific heat C vs T at 4> — 
calculated for several ratios Ao/i?Th- (b) c/Cn vs a: at = 
and T — 0.5i?Th/fcB calculated for several ratios Ao/i?Th. 
Cn = jT is the specific heat in the normal state. 



In Fig. [T{a) the electronic specific heat C vs tempera- 
ture T for a superconductor, the N region of a SNS junc- 
tion, and a normal metal are compared. In the normal 
metal case, the low-temperature specific heat, i.e., Cn = 
7T, where 7 — 27r^7\/Ffc|/3 is plotted. For the super- 
conductor, the figure displays the standard BCS result, 
where Cs follows the low-temperature exponential drop 
according to Cs(T) ~ 1.347rc(Ao/A:BT)3/2e-^«/'=BT f^^. 

T ^ Tc, and exhibits the discontinuous jump SCs = 
1.437rc at the critical temperature In the SNS 

case we have chosen a moderately-long junction with 
Aq = 10i?Th, and we set (j) = 0. It is evident that in 
the N region of a SNS Josephson junction C strongly 
departs from the normal metal behavior, and resembles 
that of a superconductor. More in particular, it shows a 
low-temperature suppression with respect to the N state, 
with vanishingly small values for T <^C Tc, and a discon- 
tinuous jump at Tc which is smaller than that in S. The 
temperature determining the onset of strong suppression 
is set by the value of the minigap. 

The dependence on position (x/L) along the N wire 
is shown in Fig. W[h), where the electronic specific heat 
per unit length c(x, T, (jj) vs T is displayed for the same 
SNS junction. This quantity is defined according to the 
expression C{T,(j)) = c(x,T, 0). First of all, 

we note that at the NS interface (i.e., for x/L = 1/2) c 
coincides with the specific heat of a BCS superconduc- 
tor. By moving toward the center of the junction, the 
discontinuous jump at T = Tc gets reduced, while the 
temperature which sets the onset of strong suppression 
becomes smaller. This can be understood by noting that 
moving away from the NS interfaces the superconducting 
correlations induced in the N region become weaker (i.e., 
the proximity effect gets suppressed), thus making the 
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specific heat to approach that in the normal state. 

The role of junction length L is shown in Fig. [UJa), 
where the SNS specific heat C is plotted as a function 
of T for several Ao/i?Th ratios eX (j) — Q. By increasing 
L, or equivalently for large Ao/i?Th ratios, C approaches 
the normal-state specific heat, which follows from the 
weakening of the proximity effect in longer SNS junc- 
tions. Moreover, the temperature which sets the onset of 
strong suppression decreases by increasing the junction 
length which stems from a reduction of the proximity- 
induced minigap in N. For instance, for AoZ-Exh = 10 at 
T = O.lTc, the specific heat turns out to be suppressed 
by about a factor of 5 as compared to that in the nor- 
mal state at the same temperature. This points out that 
moderately-long Josephson junctions are more suitable to 
maximize the effects of superconducting correlations on 
the specific heat. 

On the other hand, it is interesting to investigate the 
behavior for fixed junction length, i.e., for fixed i^xh, but 
for superconductors with different energy gap Aq. Fig- 
ure[2jb) shows the low-temperature SNS specific heat per 
unit length vs position a: at T = 0.5i?Th and ^ = for 
several values of Aq. By moving away from the middle of 
the junction {x — 0) the specific heat strongly decreases, 
approaching the value of the superconductor at the NS 
interface {x = L/2). Furthermore, by reducing the super- 
conducting gap Aq, the specific heat value in the middle 
of the N region tends to that in the normal state. For 
instance, an increase of Aq by a factor of 10, i.e, from 
10i?Th to 100i?Th which corresponds roughly to change 
from aluminum (Al) to niobium (Nb) electrodes, trans- 
lates in a reduction of C by almost a factor of 3. This can 
be understood by noting that for superconductors with 
smaller gap and fixed junction length, the proximity ef- 
fect is weakened in the middle of the Josephson junction, 
leading to an enhancement of the electronic specific heat 
value. For a given junction length, larger gap super- 
conductors are thus preferable in order to suppress the 
low-temperature specific heat. 

Figure [3l^ a) shows the SNS specific heat C versus tem- 
perature T for a junction with Aq = lOi^Th at different 
values of the phase difference 4> between the superconduc- 
tors. For a comparison, the electron specific heat in the 
N and S state are plotted as well. By increasing (f) from 
zero, the minigap in the N region DOS decreases [l3| 
therefore enhancing the low-temperature specific heat. 
At = TT (not shown) the minigap vanishes [l3| , so that 
C approaches the N state specific heat at very low tem- 
peratures (i.e., for T — > 0). 

The full phase dependence is plotted in Fig. El^b), 
which shows the specific heat in the N region of a SNS 
junction with Aq = 10i?Th at different temperatures T . 
C is a non-monotonic function of 0, and it is maximized 
at an optimal phase difference which strongly depends on 
T . In particular, the optimal value shifts toward larger 
by lowering the temperature. Notably, the figure shows 
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FIG. 3: (color online) (a) Specific heat C vs T calculated for 
Ao — lOii^Th and different values, (b) Specific heat C vs 
calculated for Aq — 10i?Th and several T values. 



that, especially at the lowest temperatures, the specific 
heat variation with can be extremely large. For in- 
stance, at T = O.lTc, the specific heat is enhanced by 
about 340% from = to ~ O.Stt, while at T = O.OSTc 
a huge variation of the order of 8800% can be obtained by 
varying the phase from to 0.97r. This demonstrates the 
sensitivity of the electronic specific heat to (f> which stems 
from the coherent nature of superconducting correlations 
induced in SNS Josephson junctions. 

We now discuss a possible experimental setup for mea- 
suring the predicted effects. A realistic measurement 
scheme is shown in Fig. |4l The basic element of the struc- 
ture consists of an ac proximity superconducting quan- 
tum interference device (SQUID), i.e., a superconducting 
loop interrupted by a normal metal region of length L. 
The electronic specific heat C in the N region can be con- 
trolled by an externally applied magnetic field B, giving 
rise to a total flux $ through the loop area. The system 
can be realized by state-of-the-art e-beam lithography 
and shadow mask evaporation of metallic thin films. An 
array consisting of a large number of identical proximity 
SQUIDs can be deposited on a silicon suspended mem- 
brane on which a thermometer and a heater have been 
previously patterned. The structure specific heat can be 
measured by ac calorimetry 20], which consists of in- 
ducing temperature oscillations of the thermally isolated 
membrane and thermometer by supplying ac power to 
the heater. 

In the following we shall provide some estimates to 
get the order of magnitude of what can be experimen- 
tally measured. For our estimate we choose a 50-nm- 
thick copper (Cu) wire, with L = 0.6 /im, width of 
200 nm, D = 0.01 rn^s''^, and density p = 8.92 x 10^ 
kgm~^. With such parameters we get a N region vol- 
ume V = 6 X 10~^^ m'^, mass m = 5.35 x 10~^^ kg, and 
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FIG. 4: (color online) A possible experimental setup consist- 
ing of a SNS ac proximity SQUID threaded by an external 
magnetic field. The loop geometry will allow to change the 
phase difference at the NS boundaries by the application of 
a magnetic field (B), thus controlling the electronic specific 
heat of the N region. An array consisting of a large number of 
identical loops (~ 10^) is patterned onto a silicon suspended 
membrane, while the heater and thermometer will allow to 
measure the structure specific heat by ac calorimetry [20l ]. 



i?Tii — 18/zeV. By choosing Al for the superconducting 
loop with Ao = 200 /ieV we get Ao/i^Th ^ 10.9. Fur- 
thermore, at r = 100 mK the copper heat capacitance 
is C^" — 2 X lO^'^Jkg^^K"^ which gives a total specific 
heat of the N island Cg" = mC§" ~ 1 x lO-^^JK-i. 
From our calculations for a SNS Josephson junction with 
Ao / -Exh = 10 and (f> — Owe have that a specific heat vari- 
ation of 1600% is predicted by changing the temperature 
from O.lTc to 0.2Tc (i.e., from about 120 to 240 mK for 
Al) [see Fig. ^a)] , and a variation of the order of 340% at 
T = O.lTc by changing (j) [see Fig. [3][b)]. Thanks to prox- 
imity effect, a change of C of 1.6 x 10~^^JK^^ by varying 
the temperature in such a range, and of 3.4 x 10~^°JK~^ 
by varying the phase is expected for the system pro- 
posed in Fig. 3) Recent developments in highly-sensitive 
calorimetry applied to mesoscopic systems allow to mea- 
sure such small signatures |2ll . [2^ . The sensitivity of 
these experiments can be estimated of the order of IO/cb 
per /.tm^ of the structured silicon membrane. A super- 
conducting loop of 1/im of diameter occupies a surface 
of approximately 4/xm^ in order to avoid any interaction 
with the neighboring loops of the array, so that a min- 
imum of ^ 5.5 X lO^^^J/K can be measured. Signal to 
noise ratios ranging from several tens to some 10^ are ex- 
pected, so that such signatures could be detected at these 
temperatures. In practice, what will be measured is the 
specific heat of the whole structure which includes the 
contribution from the membrane and transducers, and 
that related to n (typically ^ 10®) identical loops of the 
array. The loop specific heat is the sum of the N and S 
parts according to Eq. In particular, at low temper- 
ature, i.e., below about T = 0.2Tc, the N region provides 



the predominant contribution to the specific heat of the 
loop [see Fig. [3l[a)], while at higher temperature the one 
from the superconductor becomes relevant. 

In summary, we have presented a microscopic descrip- 
tion of the influence of proximity effect on the electronic 
specific heat of a mesoscopic SNS Josephson junction, 
and suggested a scheme for its experimental verification. 
We found that the specific heat can be dramatically af- 
fected by superconducting correlations. In particular, 
the strongest effects are predicted for intermediate or 
moderately-long diffusive SNS junctions, where both at 
low and high temperatures drastic deviations from the 
behavior of a normal metal are to be expected. No- 
tably, the magnitude of the electronic specific heat can 
be finely and effectively controlled by changing the phase 
difference between the superconductors. This "active" 
phase-tuning of the thermodynamic properties of meso- 
scopic conductors may pave the way for novel applica- 
tions, e.g., ultrasensitive single-photon Josephson detec- 
tors based on proximity effect ^9] . 
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